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the cluster size distribution function as generalized homogeneous function of the system size and 
two connectivity lengths. The proposed scaling theory has been verified numerically on two different 
anisotropic percolation models. 
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I. INTRODUCTION 

Percolation, a model of disorder [ij, shows a geometrical phase transition at the percola- 
tion threshold characterized by singularities of cluster related quantities similar to the critical 
phenomena or second order phase transition in thermodynamic systems |2|. The singularities 
occurred in phase transitions are generally described by power laws characterized by well 
defined critical exponents. The set of critical exponents and the scaling relations among 
them characterize the universality class of a system. Since most systems are not solvable 
analytically, numerical methods are employed in order to investigate the critical behavior. 
At the same time, the numerical results are very often limited by the finite system size. 
In a finite system, there is rounding and shifting of critical singularities depending on the 
ratio of correlation length C, to the linear dimension L of the system. In order to obtain the 
behavior of the infinite systems, the results of finite systems are generally extrapolated using 
finite-size scaling (FSS)j3||. However, there are very few FSS theory for anisotropic models 

1. FSS theory is available for isotropic percolation models such as 
5| and spiral percolation (percolation under rotational constraint) 6|. 
Anisotropic percolation clusters are generated if an external space fixed directional con- 
straint is applied on the percolation model. Such models are directed percolation (DP)[7| 
and directed spiral percolation (DSP)|8|. In these models, there are two connectivity lengths 
^11 and ^_L and they become singular at the percolation threshold with two different criti- 
cal exponents. Considering DP as a minimal stochastic Markovian process represented by 
Langevin equation, a FSS theory has been developed by Janssen et a/(9|] below the upper 
critical dimension d^ = 4: and above dc by Liibeck and Janssen |lOl|. Usually, the anisotropy 
considered in the FSS theory is either in the interaction between the constituent particles 4] 
or in the topology of the system (strip like structure) [.?]. However, in a geometrical model 
like percolation, the particles are simply occupied with a probability p in absence of any in- 
teraction between them. The models are generally defined on topologically isotropic systems 
of size L X L. For such noninteracting anisotropic models defined on topologically isotropic 
systems, a simple phenomenological FSS theory is proposed here. The proposed FSS theory 
for anisotropic percolation models is developed here considering the cluster size distribution 
function as a generalized homogeneous function [2,] of occupation probability p and ratios of 
connectivity lengths with the system size. 



II. ANISOTROPIC PERCOLATION MODELS 

There are two well known anisotropic percolation models DP[7| and DSP[8|]. In DP, a 
directional field E is present in the model. The direction of applied E field from upper left to 
the lower right corner of the lattice is considered here. As an effect of the E field, the empty 
sites to the right and to the bottom of an occupied site are only eligible for occupation. The 
eligible sites are then occupied with probability p. Accordingly, the clusters grow in the 
diagonal direction along E. In the case of DSP, a crossed rotational field B is also present 
in addition to the directional field E. In this problem, E field is applied from left to right in 
the plane of the lattice and B is applied perpendicular to E and into the plane of the lattice 
(viewed from top). Due to E field, empty site on the right of an occupied site is eligible 
for occupation whereas for B field, empty sites in the forward and clockwisely rotational 
directions are eligible for occupation. The forward direction is the direction from which the 
present site is occupied. Because of the simultaneous presence of both the E and B fields 
crossed to each other, a Hall field appears in the system perpendicular to both E and B. 
As a result, an effective directional constraint E^ff acts on the system along the left upper 
to right lower diagonal of the lattice. The clusters grow along the effective field £^eff|8|. A 
cluster is considered to be a spanning cluster, if either the lateral or the vertical extension 
of a cluster becomes equal to the dimension of the lattice L. At the percolation threshold 
Pc, a spanning cluster appears for the first time in a system. It was found as Pc ~ 0.705489 
for DP|7| and pc ~ 0.6550 for DSP[8] on the square lattice in case of site occupation. For a 
given finite system L, the cluster properties are calculated generating finite clusters. 

Typical large clusters for DP and DSP generated on 256 x 256 square lattice at their 
respective percolation thresholds are shown in FiglH It can be seen that the clusters are 
anisotropic and rarefied. In order to characterize the cluster's connectivity property, two 
length scales C,\\ and C,± are required, ^y is the extension along the elongation of the cluster 
and ^_L is the extension in the perpendicular direction of the elongation. Two connectivity 
lengths ^11 and ^± diverge with different critical exponents u^ and z/^ at the p = pc- Though 
the DSP cluster grows along an effective directional constraint E^ff, the critical properties 
of the clusters aX p = pc were found different from that of DP clusters. Accordingly DP and 



DSP belong to two different universality classes 
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III. FSS THEORY 

A system is said to be finite if tfie system size L < ^, tlie connectivity lengtli. In tfie 
case of anisotropic percolation models, there are two connectivity lengths ^y and ^± where C,\\ 
always greater than ^_|_. The finiteness of the system size then can be defined in terms of C,\\- 
According to the theory of critical phenomena, thermodynamic functions become generalized 
homogeneous functions at the critical point [2]. The cluster size distribution function Ps{p, L) 
describing the geometrical quantities here in percolation is then expected to be a generalized 
homogeneous function at the percolation threshold pc- A simple phenomenological FSS 
theory for these anisotropic percolation models is proposed here assuming Ps{p,L) as a 
generalized homogeneous function. In order to develop FSS theory, the scaling of the order 
parameter of the percolation transition Poo, probability to find a site in a spanning cluster, 
with the system size L is considered first. At the end, the scaling form will be generalized 
for arbitrary cluster related quantity Q. 

In the case of an infinite system, the order parameter Poo becomes singular at p = pc as 

Poo-^ip-Pc)'' (1) 

with a critical exponent (3. In a finite geometry with topologically symmetric dimension 
L X L, the critical singularity of Poo for anisotropic percolation clusters depends not only 
on {p — Pc) but also on the the ratios C,\\/L and ^±/L. The functional dependence of Poo on 
these parameters is then given by 

Poo(L,j>)=^[(p-p,),|,^]. (2) 

For an infinite system, at p = pc, C,\\ is infinitely large and the system properties become 
independent of the system size L. Accordingly, two parameters C,\\/L and C,±/L in Poo given 
in Eqj2]can be reduced to a single parameter ^||/^±. Thus, Poo in Eqj2]can be expressed as 

Poo = Q[{p-Pc),^\\/^±\. (3) 

It is now important to know how ^y and C,± scale with the dimension L of the finite system 
at the percolation threshold Pc- Two new scaling forms for C,\\ and C,± with L are assumed 

as, 

^11 ^ L'^ii and U ^ L'^ (4) 
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where 9\\ and 9± are two new exponents. It is now possible to define Poo in terms of tlie 
system size L as 

Poo = F[(p-p.),I^'"-'-]. (5) 

Assuming tlie order parameter Pqo as a generalized homogeneous function of {p — Pc) and 
L^ii~^^, it is possible to express Poo as 

F[X%p-p,),X''L'^^-'^] = XPoo (6) 

where a and b are arbitrary numbers and A is a parameter. The above relation is valid for 
any value of A. For A = L^^^w^^^y''^ p^ takes the form 

Poo = L^^'^^''-^F[{p-p^)L-^^h-<^-\l\ (7) 

where A = 1/h and B = a/b are two exponents to be determined. However, as L ^ oo, 
the L dependence of Pqo will vanish. Therefore, F[z] should go as z^^^ in the limit L ^ oo 
when z = {p — p^)-^"^ "^^^ • In that case, 

Po. ^ L^^'i^-'-^ {L-mr^-\p - p,)f 

^ {p-Pc)^'"'- (8) 

The order parameter exponent (3 is then given by /? = A/B. 

The size of a cluster is given by the number of occupied sites s in that cluster. For 
anisotropic clusters, there are two connectivity lengths which could be measured in terms of 
radii of gyrations R\\ and R±_ with respect to two principal axes of the cluster. The size of 
the cluster can also be calculated in terms of area defined by R\\ and P_l. It is expected that 
the cluster size should scale as s ~ P||P_|_^ at p = pc and it should go as s ~ R\\R± 
above p^ where d is the spatial dimension of the lattice and df is the fractal dimension of 
the infinite clusters generated on the same lattice. The percolation probability Poo is the 
ratio of the number of sites on the infinite cluster to the total number of sites and can be 
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given as 

,(rf/-i) 

R\\Rj_ 

Assuming P^ ^ ^_l ~ L'^^ at p^, Eq|9] leads to Poo ~ L^^i'^f-'^). Also at p = Pc, the 
functional form of Poo, given in EqjTl reduces to Poo ~ L^'-^n^^-^-'. Therefore, exponent A 



can be obtained in terms of the new exponents 9\\ and 9±_ as 

Inserting the value of A in EqJT] aX p = pc it reduces to 

P^ = L^^('^/-'^)F[0] = ei'^-"V[0] = e!''^'^'"^'"F[0] (11) 

where F[0] is a constant. For an infinite system, at p = Pci the connectivity lengths diverge 

as 

eii-b-Pci-^ii a-b-Pci"^^ (12) 

where i/y and z/^ are connectivity exponents. The following scaling relations then can easily 
be extracted as 

/3 = v^[d - df) and (3 = iy\\{d - df)e^/e\\. (13) 

The first one of these relations is the hyperscaling relation [l|] and the second one is a new 
scaling relation connecting the exponents d_[_ and 6*11. Using the above scaling relations and 
eliminating [d — df), the values of the exponents A and B can be obtained as 



A 



0±(3 _ ^1,/? 



- 0^)iy^ (^11 - e^)iy\\ 



n n 

From the expressions of A and B, it can be seen that i^w/O^ = i'±/6±. This is consistent 
with the assumption ^± ~ ^i," ^ in EqJH and the fact ^_l ~ ^m" ^ for an infinite system. 
Consequently, one can define the anisotropy exponent 

e = e\\/e^ = iy\\/iy^ (15) 

as suggested in Ref.|12|]. This is interesting to note that this equation is always valid even if 
the hyperscaling relations are not exactly satisfied. It is known that hyperscaling relations 



are not satisfied in case of DP 13j 



Since the values of A and B are now known, the finite-size scaling form of Pqo can be 
given as 

PooiL,p) = L-^^ii/'^iiF[(p-p,)L^ii/'^it] = L-'^'^^^''^F[{p-p,)L^^/''^]. (16) 



It is interesting to note that a finite size scaling relation is obtained in terms of two unknown 
exponents 6*11, 9±^ and two known exponents v\\, v±. 

Now, the finite size scaling form of the order parameter Pqo can be generalized for any 
cluster related quantity Q. Suppose that in an infinite system, as p -^ Pc the cluster related 
quantity Q scales as 

Q-b-Pcl-" (17) 

where g is a critical exponent. In a finite geometry, it is then expected that the cluster 
related quantities in general will obey the finite size scaling law given by 

g(L,p) =L^^ii/'^iiF[(p-pe)L^ii/'^ii] =L'?''^/^^F[(p-p,)L^^/^^]. (18) 

It should be mentioned here that the finite size scaling obtained here for the anisotropic 
cluster is very similar to that of isotropic clusters. In isotropic FSS, Q{L,p) = L''^'^F[(p — 
Pc)L^''^] where u is the connectivity length exponentp]. In EqlTSl 1/u is replaced by 9\\/i^\\ 
or 6±/h'±_ since it is assumed that the connectivity lengths C,\\ and C,± scale with the system 
size L with exponents 6*11 and 6±. However, it should be noticed that the FSS relation with 
respect to C,±, the transverse length, is new and nontrivial. 

The exponents 6\\ and d± are measured below for both DP and DSP clusters and the FSS 
theory developed here is applied to the cluster related quantities of both the models. 

IV. VERIFICATION OF FSS THEORY 

In order to verify the proposed FSS theory, simulations are performed on the square lattice 
of sizes L = 128 to 2048 in multiple of 2. Average has been taken over 5 x 10^ large finite clus- 
ters. First, the exponents 6\\ and 6± which describes the dependence of connectivity lengths 
with the system size L are determined. Since, clusters are grown following single cluster 
growth algorithm, the connectivity lengths are given by ^^ = 2 J2'^ R?<sPs{p, L)/ J2'^ sPs{p, L) 
and ^i = 2 X)s R±sPs{p, L)/ J^'^ sPs{p, L) where R\\ and -R_l are radii of gyration with respect 
to two principal axes of the cluster. R\\ and R_i are estimated from the eigenvalues of the 
moment of inertia tensor, a 2 x 2 matrix here. The cluster size distribution function Ps{p, L) 
is defined as Ns/Ntot where Ng is the number of s-sited clusters out of Ntot clusters generated 
on a given system size, ^y and ^_l are measured for various system sizes L at p = Pc for 
both DP and DSP clusters and plotted against the system size L in Figl2]J^a) and Figj2](6) 



respectively. The squares represent the data of DP and the circles represent that of DSP. 
There are two things to notice. First, the exponent 9\\ is found ~ 1 for both the models. 
It is expected. Because, for the given field configuration the clusters are elongated along 
the diagonal of the lattice and ^y, extension along the elongation of the cluster, then should 
be ~ \/2L for large clusters. Second, the exponent 6±^ is found different for DP and DSP: 
6^ ^ 0.64 ± 0.01 for DP and 6^ ^ 0.83 ± 0.01 DSP. Assuming e\\ = 1, one should have 
0± = ^±/^\\- Since for DP, u^ = 1.0972 ±0.0006 and u\\ = 1.7334 ±0.001 the expected value 
oi6± is ~ 0.633. Similarly for DSP, the expected value of ^^ is ~ 0.84 since i'±_ = 1.12 ±0.03 
and z/|| = 1.33 ± 0.01. The measured values are close to the expected values for both DP 
and DSP. This implies that the system properties scale with the system size L aX p = Pc 
in the similar manner as it behaves with p around pc for a given large system size. It can 
also be noticed that the magnitude of C,± is different for DP and DSP. C,i_ is larger for the 
DSP clusters in comparison to that of DP clusters whereas C,\\ of DSP clusters is slightly 
smaller than that of DP clusters. Consequently the DSP clusters are less anisotropic than 
DP clusters. 

Next, the average cluster size x is measured for different system size L aX p = pc- In 
single cluster growth algorithm, the average cluster size is defined as x = J2'sSPs{p,L), 
where Ps{p,L) is the cluster size distribution function. In an infinite system, x diverges as: 
X ~ |p — Pc\~'', 7 is a critical exponent. According to the FSS theory, it should behave as 
x(L) ~ L'^^iu'^w or x(-L) ~ L'^^^l^^ qX p = pc- In FiglHl the average cluster size x is plotted 
against the system size L for both DP and DSP. The cluster size x follows a power law with 
the system size L. The obtained slopes are 1.31 ± 0.01 for DP and 1.38 ± 0.01 for DSP. 
The expected values of the ratio of the exponents are 'y/i^w ~ 1.31 and '~f6±/i'j_ ^ 1.33 for 
DP where 7 = 2.2772 ± 0.0003. For DSP 7 = 1.85 ± 0.01 and the expected ratios of the 
exponents are: ■y/i^\\ ~ 1.39 and ■y9±/h'±_ ^ 1.37. It can be seen that the measured values 
are in agreement with that of the expected values within error bars. It confirms that the 
cluster properties follow the proposed finite size scaling theory. It is not surprising that the 
scaling theory works with ^y. However, the proposed theory works with ^_l, the transverse 
length. This is a new result. 

Finally, the FSS function form has been verified. The average cluster size is given as 
x{L,p) = L^^ii/^iiF[(p-pjL^ii/'^ii] = L^^^/''^F[{p-p^)L^^/^^] for different system size L. 
In FigHl the scaled average size x/ LP^^^ ^^^ is plotted against the scaled variable z = [p — 
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p^)L^\\/'^\\. A reasonable data collapsed is obtained for both DP and DSP. In the inset of 
FiglH the data collapse is shown by plotting x/L^^^/'^^ versus z = {p — pc)L^^^'^^. It can 
be seen that the tail of the scaling function F{z) shows a power law behavior in both the 
models with two different scaling exponents, approximately 2.23 for DP and 1.86 for DSP. 
Once again it confirms that DP and DSP follow anisotropic finite size scaling and belong 
to two different universality classes. Note that, these exponents are close to the respective 
cluster size critical exponents 7 for infinite systems as expected. 

V. CONCLUSION 

A finite size scaling theory is proposed here for anisotropic percolation models like DP and 
DSP. In this theory the cluster size distribution is assumed to be a generalized homogeneous 
function of the {p — pc) and the ratios of the connectivity lengths .^n and ^_l to the system 
size L. Reducing the functional dependence from three variables to two variables {p — Pc) 
and ^\\/^± and assuming a new scaling form of ^y and ^_l with the system size L, a finite 
size scaling form of the cluster properties are obtained. Following a single cluster growth 
Monte Carlo algorithm, clusters are generated at the percolation threshold for DP and DSP 
varying the system size L. The numerical simulation confirms the proposed scaling relations 
as well as the scaling function form. It could be considered as a simplest possible anisotropic 
finite size scaling theory for directed percolation models. 
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(a) DP 




(b) DSP 






FIG. 1: Typical large clusters of (a) directed percolation and (b) directed spiral percolation gen- 
erated at p = pc on a square lattice of size L = 256. Arrows represent the directional field E and 
the encircled dot represents the rotational field B. 
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FIG. 2: Plot of (a) ^m versus system size L and (6) ^_|_ versus L for the DP (□) and DSP (Q) clusters 
at their respective percolation thresholds. From the slopes, values of ^n and 6^ are obtained as 
6'|| = 1.01 ± 0.01 and 6^ = 0.64 ± 0.01 for DP clusters. For DSP, the values are e\\ = 1.01 ± 0.01 
and 9±_ = 0.83 it 0.01 respectively. 
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FIG. 3: Plot of average cluster size x versus system size L for the DP (□) and DSP (Q) clusters 
at pc- From the slopes, the ratio of ^Oh/vh or 7^_l/z^_l is obtained as 1.31 it 0.01 and 1.38 it 0.01 
for DP and DSP respectively. 
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FIG. 4: Plot of scaled average cluster size xiP^ L)/L'^ l'''^!! versus scaled variable \p — pc\L H'^^H for 
DP and DSP clusters for different values of L and p. The data plotted correspond to different 
system sizes of L = 128 (Q), 256 (□), 512 (O), 1024 (A) and 2048 (x). The \p-pc\ values here are 
0.01 to 0.10 in the interval of 0.01. In the inset, the data collapse is shown plotting x(p, L)/ L"'^^'^^ 
versus \p — pc\L^^'^^- 
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